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o . 

f*^ , We discuss the recent derivation of the two loop virtual QCD corrections to the W 

rsj ■ boson pair production in the quark-anti-quark-annihilation channel in the limit where 

all kinematical invariants are large compared to the mass of the W boson. In particular, 
we describe the use of the PSLQ algorithm on an example integral. 

^ ■ 1 Introduction 

^ C| The Large Hadron Collider (LHC) will be the centre of interest for particle physics phe- 

Qlh' nomcnology in the next years. Open issues that require definite answers are the verification 

of the consistency and validity of the Standard Model (SM) in the energy range of the LHC 

qj '" as well as insights into New Physics. Probably, the most important goal for the LHC is 

,Xh , the discovery of the elusive Higgs boson. Another important endeavour at the LHC is the 

precise measurement of the hadronic production of gauge boson pairs. Deviations from the 

SM predictions would indicate the presence of either anomalous couplings or new heavy 

^- ' particles which would decay into vector boson pairs [2,3]. 

In this context, W pair production via quark-anti-quark-annihilation, 

*£V qq -> W + W~ , (1) 

r~^. . is a very important process at the LHC. Firstly, it will serve as a signal process in the search 

for New Physics since it can be used to measure the vector boson trilinear couplings as 
predicted by the Standard Model (SM). Secondly, qq — > W + W~ is the dominant irreducible 
background to the promising Higgs discovery channel pp — » H — ► W*W* — » Wl'u', in the 
mass range Mniggs between 140 and 180 GeV [4]. 

Due to its importance, the study of W pair production in hadronic collisions has attracted 
a lot of attention in the literature [5-11] and is currently known at NLO. Nonetheless, if 
a theoretical estimate for the W pair production is to be compared against experimental 
measurements at the LHC, one is bound to go one order higher in the perturbative expansion, 
namely to the next-to-next-to-leading order (NNLO). This would allow, in principle, an 
accuracy at the level of 10%. High accuracy for the W pair production is also needed when 
the process is studied as background to Higgs discovery channel gg — > H. The signal in this 
case is currently known at NNLO level [12-22]. Another process that needs to be included 
in the background is the W pair production in the loop induced gluon fusion channel, 
gg — > W + W~ which contributes at 0(a 2 s ) relative to the quark-anti-quark-annihilation 
channel but is nevertheless enhanced due to the large gluon flux at the LHC [23,24]. 

Here we discuss some technical details occurring in the computation of the NNLO virtual 
corrections for W pair production [25-27]. Our methodology for obtaining the massive 
amplitude (massless fermion-boson scattering was studied in [28]) is very similar to the one 
followed in [29-31] which is, at its turn, an evolution of the methods employed in [32,33]. 
The amplitude is reduced to an expression that only contains a small number of integrals 
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Figure 1: A non-planar two-loop master integral. 



(master integrals) with the help of the Laporta algorithm [34] . Next comes the construction, 
in a fully automatised way, of the Mcllin-Barncs (MB) representations [35, 36] of all the 
master integrals by using the Misrepresentation package [37]. The representations are 
then analytically continued in the number of space-time dimensions by means of the MB 
package [38], thus revealing the full singularity structure. An asymptotic expansion in the 
mass parameter is performed by closing contours and the integrals are finally resummed, 
either with the help of XSummer [39] or the PSLQ algorithm [40]. In this paper we give an 
example of the use of the latter. 



2 The Calculation 

Consider the master integral shown in Fig. [T] The MB representation is 



/zoo ' 
J dz d (-m 2 y 2c ~ Z1 - 2 (- s )-*3 ( _ t) z 1+ * 2 +*3 ( _ u) -* 2 
-ioo j=1 

xr(z 2 )r(-z 1 -Z2- z 3 )r(z 3 )r{ Zl + z 3 + i)r(-z 4 )r(-z 2 + z 4 + z 5 + 1) 

xT(-2e + z 2 + z 3 - z 4 - z 5 - l)r(-z 5 )r(-2e - z\ — Zi — z 3 - z± — z% — 1) 

xr(-z 6 )r(zi + z 2 + z 4 + z 6 + l)T(-z 3 + z 5 + z 6 + l)r(-2e + z 4 - z 7 ) 

xT(e + Z\ + Zi + 2 5 + z 6 + 2)r(-3e — z 4 - z 4 — z 5 - z 6 - z 7 - 2)T(-z 7 ) 

xT(-e - z 4 + z 7 )T(2e + z\ + z 4 + z 5 + z 6 + z 7 + 2) 

x (r(-3e)r(-4e - zi - 2z 4 -z 5 -z 6 - 2)T(zi + z 5 + z 6 + 2) 

x T(zi + 2z 4 + z 5 + z Q + 2)r(-2e - z 7 )y x , (2) 

where, s, t, u the Mandclstam variables and m the mass of the W. 

One of the most important steps in our calculation is the computation of the master in- 
tegrals in closed analytic form after expansion in the mass. For that, one needs in most cases 
to resum integrals that are constants and which appear after reducing the dimensionality of 
the initial MB representations by applying Barnes' Lemmas. The constants are expressed 
in terms of the transcendental basis {l, C2 = tt 2 /6, C3j C4j"- 

While performing the steps described in the last paragraph of the Introduction for the 
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computation of I mas tcrj one needs to resum the following integral 

/•ci+ioo pc 2 +ico -i 

Iconst = / dzi I dz 2 —, -T (r(l - Zi)r(zi)r(l - Z 2 )T(-Z 1 - Z 2 + 1) 

J Cl -ioo Jc 2 ~ioo 4 l z 2 - L)Z 2 

T(z 2 )T(z 1 + z 2 ) (-^(- Zl ) - 4(z x - z 2 )^°H- Zl -Z2 + 1) 

+4zi^ (0) (-zi - z 2 + 2) - 4z 2 ^ (0) (-zi -z 2 + 2) + ^ (0) («j - 1))) , (3) 

where ip(°\x) is the digamma function, c\ = —1/3, c 2 = —1/3 and the integration contours 
are straight lines. 

Using a version of the double exponential algorithm [41], it takes about five minutes to 
integrate numerically Eq. [3] and get to about 30 digits. 

I™ t = 0.6265076409145894496726188621488 . (4) 

Next, using the PSLQ algorithm we fit I"™ t in terms of our basis and we finally obtain 

I __1 7tt 2 _5C(3) 

iconst- 2 + 24 2 ' { ' 

The numerical value of the resummed I cons t calculated with 32 significant digits accuracy is 
N[I const ,32] = -0.62650764091458944967261886214808 (6) 

and one can see that it agrees with the value in Eq. [4] for the first 30 digits. 

3 Conclusions and Outlook 

We have discussed some details of the computation of the two-loop NNLO QCD virtual 
corrections for the process qq — > W + W~ in the limit of small vector boson mass. Our main 
result was presented in [26]. This was a first step towards the complete evaluation of the 
virtual corrections. In a forthcoming publication, we will derive a series expansion in the 
mass and integrate the result numerically to recover the full mass dependence, similarly to 
what has been done in [42]. 

To complete the NNLO project one still needs to consider 2^3 real-virtual contributions 
and 2 — > 4 real ones. The real-virtual corrections are known from the NLO studies on 
WW + jet production in [43,44]. 
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